Polymer packaging and ejection in viral capsids: shape matters 
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We use a mesoscale simulation approach to explore the impact of different capsid geometries on the packaging 
and ejection dynamics of polymers of different flexibility. We find that both packing and ejection times are faster 
for flexible polymers. For such polymers a sphere packs more quickly and ejects more slowly than an ellipsoid. 
For semiflexible polymers, however, the case relevant to DNA, a sphere both packs and ejects more easily. We 
interpret our results by considering both the thermodynamics and the relaxational dynamics of the polymers. 
The predictions could be tested with bio-mimetic experiments with synthetic polymers inside artificial vesicles. 
Our results suggest that phages may have evolved to be roughly spherical in shape to optimise the speed of 
genome ejection, which is the first stage in infection. 

PACS numbers: 87.15.-v, 82.35.Lr 



In this paper we study the packaging and ejection of poly- 
mers of different flexibility into, and from, spherical and el- 
lipsoidal capsid shells. This is a model system for bacterio- 
phages which consist of a semiflexible polymer DNA (the 
genome) packaged into a rigid container (the phage capsid) 
|l ] . This system has recently attracted considerable theoreti- 
cal attention ||llllillllllEllE3l- Here we use a numeri- 
cal approach which was developed in fT7] where it was shown 
to reproduce the pauses during packing which have been ob- 
served experimentally |12]. Our main focus in the present 
work is on DNA ejection and on the impact of different cap- 
sid geometries on DNA packing and release. These issues had 
remained virtually unexplored via simulations until now. 

The packaged DNA is subject to strong energetic and en- 
tropic penalties because it is contained within a capsid whose 
dimensions are typically smaller than the DNA persistence 
length ~ 50 nm 

111 [3. This builds up an enormous internal 
pressure ~ tens of atmospheres which the viruses or bacte- 
riophages exploit to provide the simplest of attack strategies. 
Typically bacteriophages land on the surface of a bacterium 
and eject their genome into the host cytoplasm simply by tak- 
ing advantage of the internal pressure which pushes the DNA 
out of the phage once the capsid is opened. 

The diversity in the naturally occurring shapes of viral cap- 
sids is remarkable |14, 15]. The shells of the phage DNA 
which infects prokaryots, like E. coli or B. subtilis, are spher- 
ical or quasi-spherical, stiff shells. For example, in the 029 
phage, the DNA is fed into a 54 by 42 nm icosahedral cap- 
sid iflol . On the other hand, viruses infecting higher eukary- 
ots, which rely on a more complicated infection strategy than 
simple pressure-driven ejection, often have strikingly differ- 
ent, much more elongated shapes. An example here is the 
influenza virus which may be e.g. ~ 250 nm long and ~ 100 
nm wide 1 17]. Moreover, data on the internal volume of cap- 
sids, although sketchy, suggest that common spherical phages 
like T7, A and HK97 pack their genome at a density which is 

10 — 20 % larger than that encountered for slightly aspheri- 
cal phages like 029 or T4 (see Table 2 in Ref. |2]). Our results 
suggest a possible explanation for these observations. 



Work in this area is particularly timely because in vitro sin- 
gle molecule experiments have led to significant quantitative 
insights on the dynamics of packaging and ejection in vivo. 
Smith et al have measured the rate of packaging and the force 
of the motor for the 029 bacteriophage f 12j, while more re- 
cent experiments have characterized DNA ejection from the 
T5 and A capsids Both the ejection and packaging 

rates have been shown to vary consistently and reproducibly 
during the various stages of these processes. 

On the theoretical side, the effects of genome stiffness, ex- 
cluded volume and electrostatics on the DNA packaging pro- 
cess have been investigated by thermodynamic theories and 
simulations HH SM II Ell US). DNA ejection has also 
recently attracted a lot of attention among theorists: in par- 
ticular the roles of the buffer in DNA ejection experiments 
in vitro |8], the relation of ejection to translocation |9] and 
ratchet |21] models have been the subject of recent studies. 
Earlier work 1 10] pointed out that quasi-static analytic theo- 
ries for DNA release require an assumption for the underlying 
main mechanism leading to friction during ejection. 

In this paper we use the stochastic rotation dynamics sim- 
ulation model 1 22] to compare the way in which flexible 
and semiflexible polymers are packed into, and ejected from, 
spherical and ellipsoidal capsids of the same internal volume. 
Novel to this work are the explicit simulations of the ejection 
kinetics, which correctly capture non-equilibrium effects, and 
the focus on the impact of capsid geometry on the physics of 
DNA packing and releasing. We find that the slower relax- 
ation times of the semiflexible chains leads to slower packing 
and ejection rates. A flexible polymer is ejected more quickly 
from an ellipsoidal shell than a spherical one. However, at 
first sight surprisingly, this situation is reversed for a semi- 
flexible chain, which is ejected more quickly from a spher- 
ical shell. We argue that this is a consequence of balances 
between the thermodynamic force driving ejection and the 
ease with which the polymer can come to equilibrium within 
the confined space in the capsid. Recent advances in sin- 
gle molecule micromanipulation techniques and in DNA 
ejection imaging and analysis CI El put an experimental 
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verification of these predictions within reach. 

The polymer is a coarse-grained chain of = 100 beads 
joined by FENE springs, interacting via a potential V, 

4e^i (cr/l fi - Ti-i 1)^^ + /^i;i(ri+i - Ti) • {ri - ri-i) (1) 

where fl is the position of the i^^ bead. The first term is the re- 
pulsive part of a Lennard- Jones potential which generates ex- 
cluded volume interactions between the beads. This is in tune 
with the experiments where repulsive interactions dominate. 
The potential parameters used were e = /c^T and a =2.5nm. 
K, in the second term in eq. {I} is a bending rigidity which 
sets the persistence length / ~ ncr/kBT. Here we set / = 
for a flexible polymer and / = lOa for a semiflexible poly- 
mer. (We use 10 cr to compromise between reaching typical 
genomic stiffness - 20 a under physiological conditions LI2I] 
- and feasible length and time scales in the simulations.) The 
updating of the beads' positions and velocities is performed 
using the velocity-Verlet molecular dynamics algorithm. 
The capsid shapes, illustrated in Figure[I] are described by 

/ ^ 1 _ [(^/a)2 + (y/5)2 + (^/^)2] ^ Q 

We choose a = h = c = 3.02cr to model a sphere and a = 
h = 2.6(7 and c = 4.07cr for an ellipsoid. Each capsid is 
modeled as a hard shell with a hole that permits the entrance 
of one bead at a time. A repulsive force ksT/ (crf^) is applied 
to any bead which is at a point for which | / |< /o, where 
/o = 0.2 is a threshold. Our choices of a, b, c, and of /o lead 
to the same volume available to the chain for both shapes, 
which corresponds to a packing fraction of 0.4, comparable 
with previous numerical work and typical phage densities |3, 
|4]. Qualitatively similar results to the ones reported below 
have been found with N = 80 and = 120 with these capsid 
geometries, and with N = 200 and a packing fraction of 0.4. 

The motor that feeds the polymer into the capsid is, in re- 
ality, extremely complex jjlj]. Here we use a simple model 
aimed at capturing the basic physics. Essentially the motor 
has to (1) capture a bead and (2) feed it into the capsid. This 
is accomplished by requiring the motor to apply a radial force 
(of magnitude 5 ksT / a) if the bead enters a cylinder of ra- 
dius 0.7a and length a with origin at the capsid entrance. 
The details of this mechanism do not affect our results. Once 
captured, the bead is packed by applying a constant force to- 
wards the centre of the capsid. Our simulations allow us to 
identify the minimum motor force which is needed to achieve 
full packing. To estimate this, we ran a set of packing sim- 
ulations at different motor force, and picked the lowest value 
of the force which still, on average, packed the whole chain. 
Flexible and semiflexible polymers are respectively found to 
require a minimum motor force of 16 and 20 ksT / cr to be 
packed into a sphere, and of 18 and 26 /c^T/cr to be packed 
into an ellipsoid. In general the difference between the forces 
corresponding to the sphere and the ellipsoid increases with 
packing fraction. 

The polymer is coupled to a coarse-grained solvent model, 
stochastic rotation dynamics. This acts as a hydrodynamic 



(A) 



(B) 




FIG. 1: Schematic representation of the simulation. A semiflexible 
polymer is first packed into and then ejected from a rigid capsid. The 
first row (A) shows the capsid and the dangling tail of the polymer 
close to the end of a packing run. The second row (B) shows the 
configuration of the polymer chain inside the capsid. Column (I) 
refers to a spherical capsid, and column (II) to an ellipsoidal one. 



thermostat allowing momentum transfer between beads and 
allowing flows to be set up in the surrounding fluid as a con- 
sequence of the bead motion. The solvent has a viscosity ~ 5 
cP, comparable to that of cytosol. The capsid is permeable to 
the solvent, which is the physical situation for phage capsids. 
(We measure force and time in simulation units in Figures 2- 
4. One time and force simulation unit corresponds to 3 ns and 
1.64 pN respectively.) 

The polymer is initially configured randomly except for the 
requirement that the first bead lies within the capsid and the 
rest outside. The polymer is then equilibrated in this position 
before opening the bead entrance and applying the feeding 
force. A single bead is left out to initiate ejection once the 
motor force is set to zero. This is done after leaving time for 
the polymer to equilibrate within the capsid. 

Our simulations allow us to compare packing and ejection, 
a flexible and semiflexible polymer, and a spherical and ellip- 
soidal capsid. Figure [1] illustrates typical packed configura- 
tions for the semiflexible polymer - the polymers are ordered 
in spool-like domains (although not in an ideal inverse spool) 
as predicted theoretically |3, 6]. 

Figure |2| shows the number of packed beads as a function 
of time for both packing and ejection for the different chain 
flexibilities and capsid shapes. The motor force during pack- 
ing was 26 /c^T/cr - the minimum force to pack the semi- 
flexible polymer into the ellipsoid, the case for which pack- 
ing is hardest. The most immediately striking feature is that 
packing and ejection times are considerably faster for flexi- 
ble polymers. This is because relaxation times increase as the 
polymer becomes stiff er and e.g. once one bead has escaped it 
takes longer for the semiflexible chain to readjust itself so that 
a subsequent bead is in a position to escape. 

Perhaps more surprisingly for flexible polymers the sphere 
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FIG. 2: Number of packed beads versus time during (a) packaging 
and (b) ejection for flexible and semiflexible polymers comparing a 
spherical and an ellipsoidal capsid. 



ejects more slowly than the ellipsoid whereas for semiflexible 
polymers the reverse is true. The packing, on the other hand, 
proceeds more easily (it requires a smaller minimum force) 
and more quickly (at equal packaging force it takes less time) 
into the sphere whatever the flexibility of the polymer. 

These results can be explained by considering a combi- 
nation of entropic and dynamic arguments. Both flexible 
and semiflexible polymers lose more entropy when they are 
packed in an ellipsoid than when they are packed in a sphere. 
Therefore, based on entropic arguments alone, we would ex- 
pect the polymers to pack more easily in a sphere and to be 
ejected more easily from an ellipsoid. This is the case for the 
flexible polymers |23]. 

For the semiflexible polymers, however, the sphere is faster 
both in packing and in ejection. We believe this to be pri- 
marily a consequence of two effects. Firstly the beads of the 
semiflexible polymer suffer more from the constraint of be- 
ing in the ellipsoid when they try to rearrange themselves as 
beads are ejected. Secondly the bending energy lost in pack- 
ing is larger for the case of a sphere. These effects win over 
the entropic arguments and for the semiflexible polymer there 
is quite a pronounced advantage in ejection time for the sphere 
12411 . The difference in ejection times between the sphere and 
the ellipsoid increases with the aspect ratio of the ellipsoid. 

To investigate these non-equilibrium effects further we plot. 



FIG. 3: Force opposing the motor (in units of ksT/ a) during pack- 
aging and ejection for a semiflexible polymer in (top) a spherical and 
(bottom) an ellipsoidal capsid. 

in Figure [51 the force opposing the motor during the packag- 
ing and ejection of a semiflexible polymer for a sphere and an 
ellipsoid. We define such a force as the one felt by the bead 
inside the capsid which is closest to the motor at a given time. 
It includes the force due to local bending at the capsid en- 
trance, the elastic force due to the springs acting on the bead 
under observation, and the overall Lennard- Jones repulsion of 
the other beads, both in the capsid and in the tail outside (i.e. 
we do not include the capsid contribution). 

For both the capsid shapes there is hysteresis, i.e. the force 
during the packing is larger than the one felt during ejection. 
This shows that a significant portion of the resistance the mo- 
tor has to overcome during packing is due to dynamic dissipa- 
tive effects. The hysteresis is larger for the ellipsoidal capsid, 
supporting the assertion that chain rearrangements are more 
difficult in this case presumably due to the narrower space 
close to the capsid tip. That dynamic effects are important 
in our simulations can also be appreciated by noting that a set 
of simulations considering ejection from an ellipsoidal capsid 
with a hole on the long side yielded an ejection time compa- 
rable to that found from the sphere. 

Typically experiments report packing rates as a function of 
the number of packed beads 1 12] and we therefore present 
similar data for ejection in Figure [4^. The ejection rate de- 
creases as the number of packed beads decreases for the flex- 
ible chain. This is because the force driving the ejection is 
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FIG. 4: (a) Ejection rate as a function of the number of packed beads 
for flexible and semiflexible polymers comparing a spherical and an 
ellipsoidal capsid. (b) Data from one individual simulation of the 
ejection of a semiflexible polymer from a sphere. Arrows indicate 
pauses (corresponding to rearrangements of the polymer inside the 
capsid). 



the entropic penalty of confinement which decreases with de- 
creasing packing fraction. For the semiflexible chain, the rate 
decreases appreciably less quickly suggesting that entropic 
considerations are less dominant. One can speculate that the 
decreasing entropic force is offset by easier rearrangements 
within the capsid as it empties. 

The data in Figures |2j) and|4t are averaged over many runs 
and thus the curves appear continuous. However, individual 
runs indicate that, just as for packing, there are pauses in ejec- 
tion as the polymer rearranges itself within the capsid (Figure 

Et). 

In summary, we have compared the packing and ejection 
of flexible and semiflexible polymers within phage capsids of 
variable geometry. The behaviour is influenced not only by 
thermodynamic considerations but also by the relaxation time 
of the polymers as they try to rearrange themselves within the 
confined capsid geometry. For semiflexible polymers we find 
that a spherical shape leads to fastest ejection. One might 
speculate that this is one of the reasons why most phages 
have evolved to be roughly spherical in shape. We note how- 
ever that other phages, not relying on pressure to eject their 
genome |25], are still spherical, thus suggesting that the me- 
chanics and energetics of the proteins making up a viral capsid 
are equally important in ultimately determining shape. We 
also found that non-equilibrium effects hamper packing into 



an ellipsoidal phage, which may explain why aspherical cap- 
sids contain less DNA. Our simulations suggest a series of sin- 
gle molecule bio-mimetic experiments in which the dynam- 
ics of polymers of variable flexibility undergoing packing- 
ejection cycles into and out of vescicles of controlled shape 
are studied. 
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